BRCA1 frameshift variants leading to extended incorrect protein C termini.

Carriers of BRCA1 germline pathogenic variants are at substantially higher risk of developing breast and ovarian cancer than the general population. Accurate identification of at-risk individuals is crucial for risk stratification and the implementation of targeted preventive and therapeutic interventions. Despite significant progress in variant classification efforts, a sizable portion of reported BRCA1 variants remain as variants of uncertain clinical significance (VUSs). Variants leading to premature protein termination and loss of essential functional domains are typically classified as pathogenic. However, the impact of frameshift variants that result in an extended incorrect terminus is not clear. Using validated functional assays, we conducted a systematic functional assessment of 17 previously reported BRCA1 extended incorrect terminus variants (EITs) and concluded that 16 constitute loss-of-function variants. This suggests that most EITs are likely to be pathogenic. However, one variant, c.5578dup, displayed a protein expression level, affinity to known binding partners, and activity in transcription and homologous recombination assays comparable to the wild-type BRCA1 protein. Twenty-three additional carriers of c.5578dup were identified at a US clinical diagnostic lab and assessed using a family history likelihood model providing, in combination with the functional data, a likely benign interpretation. These results, consistent with family history data in the current study and available data from ClinVar, indicate that most, but not all, BRCA1 variants leading to an extended incorrect terminus constitute loss-of-function variants and underscore the need for comprehensive assessment of individual variants.

NERDSS: A Nonequilibrium Simulator for Multibody Self-Assembly at the Cellular Scale.

Currently, a significant barrier to building predictive models of cellular self-assembly processes is that molecular models cannot capture minutes-long dynamics that couple distinct components with active processes, whereas reaction-diffusion models cannot capture structures of molecular assembly. Here, we introduce the nonequilibrium reaction-diffusion self-assembly simulator (NERDSS), which addresses this spatiotemporal resolution gap. NERDSS integrates efficient reaction-diffusion algorithms into generalized software that operates on user-defined molecules through diffusion, binding and orientation, unbinding, chemical transformations, and spatial localization. By connecting the fast processes of binding with the slow timescales of large-scale assembly, NERDSS integrates molecular resolution with reversible formation of ordered, multisubunit complexes. NERDSS encodes models using rule-based formatting languages to facilitate model portability, usability, and reproducibility. Applying NERDSS to steps in clathrin-mediated endocytosis, we design multicomponent systems that can form lattices in solution or on the membrane, and we predict how stochastic but localized dephosphorylation of membrane lipids can drive lattice disassembly. The NERDSS simulations reveal the spatial constraints on lattice growth and the role of membrane localization and cooperativity in nucleating assembly. By modeling viral lattice assembly and recapitulating oscillations in protein expression levels for a circadian clock model, we illustrate the adaptability of NERDSS. NERDSS simulates user-defined assembly models that were previously inaccessible to existing software tools, with broad applications to predicting self-assembly in vivo and designing high-yield assemblies in vitro.

Structurally Linked Dynamics in Lactate Dehydrogenases of Evolutionarily Distinct Species.

We present new findings about how primary and secondary structure affects the role of fast protein motions in the reaction coordinates of enzymatic reactions. Using transition path sampling and committor distribution analysis, we examined the difference in the role of these fast protein motions in the reaction coordinate of lactate dehydrogenases (LDHs) of Apicomplexa organisms Plasmodium falciparum and Cryptosporidium parvum. Having evolved separately from a common malate dehydrogenase ancestor, the two enzymes exhibit several important structural differences, notably a five-amino acid insertion in the active site loop of P. falciparum LDH. We find that these active site differences between the two organisms' LDHs likely cause a decrease in the contribution of the previously determined LDH rate-promoting vibration to the reaction coordinate of P. falciparum LDH compared to that of C. parvum LDH, specifically in the coupling of the rate-promoting vibration and the hydride transfer. This effect, while subtle, directly shows how changes in structure near the active site of LDH alter catalytically important motions. Insights provided by studying these alterations would prove to be useful in identifying LDH inhibitors that specifically target the isozymes of these parasitic organisms.

Enzymatic Kinetic Isotope Effects from First-Principles Path Sampling Calculations.

In this study, we develop and test a method to determine the rate of particle transfer and kinetic isotope effects in enzymatic reactions, specifically yeast alcohol dehydrogenase (YADH), from first-principles. Transition path sampling (TPS) and normal mode centroid dynamics (CMD) are used to simulate these enzymatic reactions without knowledge of their reaction coordinates and with the inclusion of quantum effects, such as zero-point energy and tunneling, on the transferring particle. Though previous studies have used TPS to calculate reaction rate constants in various model and real systems, it has not been applied to a system as large as YADH. The calculated primary H/D kinetic isotope effect agrees with previously reported experimental results, within experimental error. The kinetic isotope effects calculated with this method correspond to the kinetic isotope effect of the transfer event itself. The results reported here show that the kinetic isotope effects calculated from first-principles, purely for barrier passage, can be used to predict experimental kinetic isotope effects in enzymatic systems.